HU-EP- 10/64 
SFB/CPP-10-114 



Flavor symmetry breaking in lattice QCD with a mixed action 



Oliver Bar, a Maarten Golterman 6 and Yigal Shamir c 

a Department of Physics 
Humboldt University, Berlin, Germany 

b Department of Physics and Astronomy 
San Francisco State University, San Francisco, CA 94132, USA 

c Raymond and Beverly Sackler School of Physics and Astronomy 
Tel-Aviv University, Ramat Aviv, 69978 Israel 



ABSTRACT 

We study the phase structure of mixed-action QCD with two Wilson sea quarks and any 
number of chiral valence quarks (and ghosts), starting from the chiral lagrangian. A priori 
the effective theory allows for a rich phase structure, including a phase with a condensate 
made of sea and valence quarks. In such a phase, mass eigenstates would become admixtures 
of sea and valence fields, and pure-sea correlation functions would depend on the parameters 
of the valence sector, in contradiction with the actual setup of mixed-action simulations. 
Using that the spectrum of the chiral Dirac operator has a gap for nonzero quark mass we 
prove that spontaneous symmetry breaking of the flavor symmetries can only occur within 
the sea sector. This rules out a mixed condensate, and implies restrictions on the low-energy 
constants of the effective theory. 
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I. INTRODUCTION 



Dynamical lattice simulations with fermions that preserve chiral symmetry [lH3[ are ex- 
tremely time consuming. The numerical cost typically exceeds simulations with Wilson or 
staggered quarks by one or two orders of magnitude For this reason so-called mixed- 
action simulations have been proposed, referring to a setup with either Wilson or staggered 
sea quarks, and domain-wall or overlap valence quarks; we will refer to such valence quarks 
collectively as "chiral" quarks. Even though mixed-action theories are not unitary, they are 
widely believed to have the correct continuum limit. A key advantage is that the valence 
sector preserves chiral symmetry except for soft breaking by mass terms. This is particularly 
beneficial for the computation of weak matrix elements. 

Quite a few mixed-action simulations with staggered sea quarks have already been per- 
formed 0. All these simulations used the configurations generated by the MILC collabo- 
ration with Asqtad-improved staggered fermions [6[. Exploratory simulations with twisted- 
mass Wilson fermions [7] or clover fermions [§[ in the sea sector have also been reported, 
and more are expected in the near future. 

Mixed-action theories can be studied in chiral perturbation theory (ChPT) at nonzero 
lattice spacing 0-11]. In particular, the dominant source for unitarity violations can be 



studied analytically. The scalar correlator, for example, is a sensitive probe for unitarity 



violations |12|-|14j]. It has been shown that the numerical data for the scalar correlator 
agrees quite well with the predictions of ChPT (lol- 17|. lending support to the validity of 
mixed-action ChPT in describing lattice data. 

In this paper we study the phase structure of mixed-action theories with two Wilson- 
like sea quarks, and any number of chiral valence quarks (and ghosts). We consider the 
"Aoki" or "large cutoff effects (LCE)" regime, where quark masses are of order a 2 , the 
lattice spacing squared. The phase diagram for theories with Wilson fermions in the sea or 
the valence sector has been studied by various authors , and an interesting nontrivial 

phase structure has been found. In the two-flavor case, depending on the sign of a low-energy 
constant (LEC) in the chiral lagrangian [l9j], there exists either a first order phase transition 
with a nonvanishing minimal pion mass or an Aoki phase (23j. The latter is characterized 
by the spontaneous breaking of isospin and parity symmetries, with two of the three pions 
turning massless. We expect these scenarios to be present in the mixed-action theory as 
well, but a priori the phase structure might be more complicated. 

To make the discussion more concrete, we begin with the effective chiral theory for the 
case of two valence quarks. We find that, indeed, the potential of the effective theory 
allows for a richer phase structure. In particular, depending on the sign of a certain linear 
combination of LECs there is a "mixed" phase, characterized by a mixed condensate built 
out of and a valence quark. 

A mixed phase of the effective chiral theory immediately raises a paradox: The mixed 
condensate spontaneously breaks the separate sea and valence flavor symmetries to the 
diagonal subgroup. As a result, the mass eigenstates are admixtures of sea and valence 
fields. The masses themselves depend on both the sea and the valence quark masses. This 
means, for example, that the two-point function of pure-sea pion fields would become a 
superposition of exponentials, all of which depend on the valence quark mass. This putative 
situation is clearly inconsistent with the very setup of mixed-action theories. Indeed, in any 



As in the Aoki phase, this condensate breaks parity. 
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numerical simulation, by construction the valence and ghost determinants exactly cancel out, 
and pure-sea physics cannot possibly depend on the presence of valence and ghost sectors. 

In order to resolve this conundrum we turn to the underlying theory, mixed-action QCD. 
The chiral Dirac operator of the valence sector has a gap for nonzero valence mass. Following 
the arguments of Ref . [24j , we prove that none of the flavor symmetries of the valence-ghost 
sector can be broken spontaneously, for any number of valence quarks. Since a mixed 
condensate would necessarily break the valence-ghost flavor symmetry group, this rules out 
a mixed condensate. 

While showing that all is well in the underlying theory, this state of affairs calls into 
question the reliability of the effective chiral theory, since the latter appears to allow for a 
mixed condensate. Specializing once more to the case of two valences quarks, we derive a 
mass inequality in the underlying theory that constrains the mass of mixed charged pions 
to be not below the smaller of pure-sea and pure- valence charged pion masses. 2 In the 
effective theory, this mass inequality implies an inequality on a linear combination of LECs. 
The latter excludes the region in the phase diagram in which a mixed phase would occur, 
thereby preventing the effective theory from making predictions that are inconsistent with 
the underlying theory. 

In conclusion, the constrained effective theory gives a consistent description of the possible 
phase diagram of mixed-action QCD. 

In Sec. [Til we introduce the effective potential for mixed-action QCD for the case of two 
chiral valence quarks, and list its symmetries. In Sec. HIT] we study patterns of spontaneous 
symmetry breaking at the level of the effective theory, focusing on the mixed condensate. 
Some technical details are relegated to the appendices. The main results are derived in 
Sec. [TV] and Sec. [V] offers our concluding remarks. 



II. THE EFFECTIVE POTENTIAL FOR MIXED- ACTION QCD 

The chiral effective lagrangian for the mixed-action theory with Wilson sea and chiral 
valence quarks has been constructed in Refs. 0, [I3]. 3 It is written in terms of the nonlinear 
field 

/ p.xn I %<h I tli \ 

(2.1) 

Specializing to the case of two sea and two valence quarks, is a four-by-four hermitian 
matrix, a two-by-two hermitian matrix, while ZJ and u are four-by-two and two-by-four 
matrices, respectively, with Grassmann-valued entries. 4 We will label the rows and columns 
of £ as u s , d s , u v , d v , u v and d v , where u stands for up, d stands for down, and s, v stand 
for sea, valence, respectively, while the tilde indicates ghost quarks. 




2 For a review of mass inequalities in standard QCD, see Ref. 25 1 

3 For an introduction, see for example Ref. [26j. 
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The fact that the two-by-two matrix exp(2(j>/ f) is not unitary but hermitian follows from a detailed study 



of the symmetries in the ghost sector 22|, [27|, [28 1 
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To the order we are working here the potential is 
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£2 £2 

V = -— str (mE" 1 + Em) - J — a str (PjE" 1 + SPj) 

-a 2 W M str (T3ST3S- 1 ) - a 2 W^ str (P^E^PsE" 1 + SP S EP S ) 

-a 2 W^ (str(P 5 S- 1 + SP s )) 2 -a 2 ^ (str (P^E" 1 - EP 5 )) 2 , (2.2) 

where str denotes the supertrace. Note that we need E -1 instead of because the ghost 
part of this nonlinear field is not unitary. Furthermore, we have that sdet(E) = 1, with sdet 
the super-determinant. The parameters m and a are proportional to the quark mass matrix 
and the lattice spacing, respectively: 

m = 2P M, a = 2W a , (2.3) 

where M is the real diagonal mass matrix. The parameters / and Bq are the familiar LECs 
of continuum ChPT, while W , Wm, W! are additional LECs associated with a nonvanishing 
lattice spacing. P$ is the projector on the sea quark sector, with P$ = (l + T 3 )/2. Explicitly, 
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M = m v l , P s = , T 3 = -1 , (2.4) 



where is the two-by-two null matrix and 1 is the two-by-two unit matrix. m s is the sea- 
quark mass, m v is the valence-quark mass, and m g = \m v \ is the ghost-quark mass, which we 
always choose equal in magnitude to the valence-quark mass. 5 Since the valence determinant 
does not depend on the sign of m v , the ghost determinant cancels the valence determinant 
exactly for this choice. We assume isospin symmetry in both the sea and the valence sector. 

In the following, we will absorb the contribution of order a in the potential into the 
definition of the sea quark mass, which amounts to shifting m s — > m s + clWq/Bq. This is not 
only convenient, but also justified, since the order-a shift contributes only to the additive 
mass renormalization of the sea-quark mass. 

In order to simplify the notation we rescale the potential by V — > Vf 2 B /4, and we write 

V = -str (ME" 1 + EM) - c x str (T 3 ET 3 E _1 ) - c 2 str (PjE^PsE- 1 + EP S EP S ) 

-c 3 (str (PsE- 1 + EP5) ) 2 - c 4 (str (P^ 1 - EP 5 )) 2 , (2.5) 

where the new coefficients q, i = 1,2,3,4, are proportional to the coefficients Wm and W(, 
i = 6,7, 8. Note that we also absorbed the factor a 2 into the q, so that these LECs are now 
of order a 2 . 

In writing Eq. (12. 5p . we implicitly assume that all terms in V are of the same order of 
magnitude and equally important for the potential. In other words, we assume that we are in 
the Aoki-regime [2l| with m v and the (shifted) mass m s both of order a 2 in magnitude. Terms 
not shown in Eq. (I2.5P are at least of order a 3 , ma ~ a 3 or m 2 ~ a 4 , with m ~ m v ~ m s . 



5 Note that m g has to be taken positive in order for the QCD path integral in the ghost sector to be 
convergent. Consistently, the sign of the ghost quark mass cannot be changed by a chiral symmetry 
transformation 1221. 
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Mixed-action QCD, and thus the effective lagrangian and its potential, are invariant under 
independent flavor rotations in the sea and in the valence sector. The symmetry group G 
has the structure [sj 

G = G S ca ® G va i , G sca = U(2) v , G val = U(2|2) y . (2.6) 

There is no symmetry connecting the sea and the valence sectors, a consequence of the dif- 
ferent fermion formulations used in each sector in the underlying lattice theory. In the limit 
of vanishing lattice spacing this group enlarges to the U(4|2) y symmetry of the partially- 
quenched continuum theory [2^] . For a vanishing valence-quark mass the symmetry group is 
larger because of the exact chiral symmetry in the valence sector. For a detailed description 



of the full chiral symmetry group in the valence and ghost sectors, see Refs. [22|, |28 



III. PATTERNS OF SPONTANEOUS SYMMETRY BREAKING 

To begin our analysis, it is useful to explore the possible patterns of spontaneous sym- 
metry breaking. We will first establish that nothing interesting happens in the ghost sector, 
and then discuss possible symmetry breaking patterns in the sea and valence sectors, in 
order to provide a context for the results that will follow in the subsequent section. 

A. Ghost sector 



The most general ground state has the form 



id 



(3.1) 



where S g is a four-by-four matrix and £ g is a two-by-two matrix. A question that im- 
mediately arises is whether the Grassmann parts u and 00 can acquire any nonzero vac- 
uum expectation values. Intuitively, one would think that this cannot happen, because a 
Grassmann- valued scale does not exist. In App. lAlwe will show that, indeed, co = to = 0. 

The constraint sdet(E vac ) = 1 implies that both submatrices S q and £ g are regular and 
that det(E q ) = det(£ g ). To the order we are working, the block- diagonal form of the vacuum 
expectation value S vac implies that the ghost and quark sectors decouple, and the potential 
is a sum of two terms, V = V q + V g . The two-by-two matrix £ g is hermitian and positive, 
and can be diagonalized by an isospin transformation in the ghost sector. This means that 
E g can be written as 

v -t 

, p 2 

with P12 > 0. In terms of these two eigenvalues, the ghost-sector part of the potential then 
becomes equal to 

V g = \m v \{ Pl + p^ 1 +p 2 + P2 1 ) + 2 Cl , (3.3) 

which is minimized by pi = p 2 = 1, or, equivalently, S g = 1. 

Thus, the effective theory predicts that isospin in the ghost sector is unbroken. 6 This is 
in agreement with the general result we will establish in the underlying theory in Sec. II V Al 
As we will see next, the situation in the sea and valence sectors is more subtle. 



^ = ( n ° ) > (3-2) 



This result generalizes to any number of flavors in the ghost sector. 
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B. Scenarios in the sea-valence sector 

We begin by introducing a convenient parametrization of the most general vacuum state 
in the sea- valence sector. From the previous subsection we have that det(E g ) = 1, hence 
det(S q ) = 1 as well; E q is thus an element of SU(4). For future use, we subdivide E q into 
blocks of two-by-two matrices, 

E q = ( f> ) . (3.4) 

Furthermore, if £ OT = p v l is proportional to the unit matrix (with p v a complex number), 
E q can be written in the form 

£ q = e** ( f^ D ) , D = exp(^r 3 /2) , £ = V^JpI 2 , (3-5) 

with = mod tt/2, \p v \ < 1, and in which r 3 is the third Pauli matrix. The proof is given 
in App. EE 

Symmetries which might be spontaneously broken by the ground state are the flavor sym- 
metry G, given in Eq. (I2.6p . and the discrete symmetries parity (P) and charge conjugation 
(C), under which the field £ transforms according to 

E XT 1 , E E T . (3.6) 

The vacuum state (13. 5p preserves charge conjugation. 7 Parity is broken by the Aoki con- 
densate, which corresponds to 9 ^ 0. The mixed condensate, which corresponds to p^ ^ 0, 
or, equivalently, \p v \ < 1, breaks parity too. 

We now turn to a more detailed discussion of possible phases. The trivial vacuum is 
parametrized by E ss = E^ = 1 and H sv = T, vs = 0. A first nontrivial example is provided 
by a vacuum expectation value of the form 

^=( e 7i)' < 37 > 

with a nonvanishing isospin condensate in the sea sector (this corresponds to p v = 1 and 
= in Eq. (13. 5p ). This vacuum state corresponds to the Aoki phase, with spontaneous 
breaking of parity and flavor [23]. The nonsinglet flavor group SU(2) sea breaks down to 
U(l). Associated with this breaking are two massless Goldstone bosons, the charged pions 
iif s . The possibility of the Aoki phase in the sea sector is expected, of course. It has been 
shown in Ref. [l9| that the existence of the Aoki phase is one of two possible scenarios for 
unquenched lattice QCD with two flavors of Wilson fermions, which is precisely the sea 
sector of the mixed-action theory we study here. 

A new phase, unique to the mixed-action theory, would be a phase in which H sv and E„ s 
are nonzero, which corresponds to a condensate mixing sea and valence quarks. One way to 
explore whether such a phase might occur is to consider the meson masses obtained when 



7 For other orientations of the vacuum, charge conjugation has a more complicated form that involves a 
flavor rotation. 
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the potential is expanded around the trivial vacuum, S q = 1. Assuming m S} m v > for 
the rest of this subsection, and expanding the potential to quadratic order in 0, we find the 
following tree- level masses for mesons made out of two sea quarks (M ss ), two valence quarks 
(M vv ), and one sea and one valence quark (M sv ): 



Ml = B (2m s + 8c 2 + 16c 3 ) = B {2m s + 8c' 2 



2) > 



M 2 VV = 2B m v , (3.8) 
Ml = B (m s + m v + 4 Cl + 2c 2 + 8c 3 ) = B (m s + m v + 4 C ; + 2c' 2 ) , 



where 



c[ = ci + c 3 , c' 2 = c 2 + 2c 3 . (3.9) 

We note that c 4 does not contribute; it only contributes to the sea r\ mass through a term 
quadratic in r\ ss . 

First, the valence meson mass vanishes when m v = 0, consistent with the fact that the 
valence sector has exact chiral symmetry. Second, from the expression for M ss , we see 
that spontaneous symmetry breaking should take place in the sea sector when c' 2 < 0, and 
m s < 4|c' 2 |, because this would drive Ml to a negative value. This is in agreement with the 



ChPT argument of Ref. [19[ for the existence of an Aoki phase. 

A new type of phase is suggested by the third equation of Eq. (I3.8p . For 2c[ + c' 2 < 0, 
Ml becomes negative when m s + m v is small enough. The negative curvature at the origin 
of field space indicates that a mixed condensate develops, alongside with mixed Goldstone 
bosons. 

Let us explore this possibility, assuming that c' 2 > 0, but 2c[ + c' 2 < 0. Since c' 2 > there 
is no Aoki condensate in the sea sector. Furthermore, let us assume that the vacuum takes 
the form 



with real p. (This corresponds to choosing all phases equal to zero in Eq. (13.51) .) The 
effective potential then reduces to 

V = -V - ci = -K + m s )p - {2c[ + c' 2 )p 2 . (3.11) 

If indeed 2c[ + c 2 < 0, we find a minimum at 

m v + m s 

p= , , , , , 3.12 
H 2\2c[ + c' 2 \ ' V ; 

provided that the right-hand side is smaller than one (otherwise, the minimum is at p = 1). 

Before we turn to the mass spectrum of this mixed phase, let us recall the familiar 
situation on the trivial vacuum. The nonlinear field is parametrized as E = exp((2i//)0), 
where the pseudoscalar field <p is expanded as 



V2(j) 



( Vss + 7T S ° S 


v / 2tt+ 


Vvs + < s 


v / 2tt+ 


\ 


V27T- 


Vss ~ K° ss 




Vvs - 




Vsv + 7T^ 


V / 27T+ 


Vw ~r Tt vv 


v / 2vr+ 




V V2 7T- 


Vsv ~ T^sv 


\/2vr- 


Vvv ~ T^ vv 


J 



(3.13) 
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in self-explanatory notation. 8 Disregarding the ghost sector, the flavor symmetry consists of 
the direct product SU(2) sca <g) SU(2) va i <g> U(l) sea -vai- 9 Pure-sea fields reside in the upper-left 
two-by-two block. They transform only under SU(2) sea and not under SU(2) va i. 

The mixed phase is realized by an expectation value for neutral, mixed meson fields, 

(Vvs) = (Vsv) = (f/V2)C , cosC = p ■ (3.14) 

This breaks the flavor symmetry spontaneously to the diagonal subgroup SU(2) sca=val . As- 
sociated with the four broken generators are the mixed Goldstone pions 

?]vs Vsvi ^vs ^sv ^vs ^s-ui ""si) — ^vs- (3.15) 

The nonzero-mass eigenstates are admixtures of fields with different transformation proper- 
ties under the separate sea and valence flavor groups. The nonvanishing masses depend on 
both m s and m v . The dependence is explicit, as well as implicit via p in Eq. (13.121) . 

Now we are facing a paradox: If, for example, we calculate the two-point function 
( 7r i(0) 7r ^( a; ))) we fi n d that it is a superposition of exponentials coming from the various 
nonzero mass eigenstates of the mixed phase. These masses all depend on m v , 10 as does the 
correlation function itself. But this cannot possibly be correct, because, by the very setup 
of mixed-action theories, the sea sector does not depend on the valence part of the action 
at all. 

In fact, this observation is a little more subtle than it appears, because spontaneous 
symmetry breaking takes place in the thermodynamical limit, whereas numerical simulations 
are always done in finite volume. 

In ChPT terminology, the analysis of the potential we have just carried out corresponds 
to being in the p- regime for the Goldstone pions of Eq. (13.151) . In order to stay in the p- 
regime in finite volume, one would have to turn on "seeds" for the given symmetry-breaking 
pattern. These would take the form of mixed mass terms that couple the sea and valence 
quarks. Were such mass terms to be introduced at the quark level in the underlying theory, 
the valence and ghost determinants would no longer cancel each other. The separation into 
sea and valence sectors would no longer apply, and there would be nothing a-priori wrong 
with finding that properties of what used to be the sea sector now depend on parameters of 
what used to be the valence sector. 

Numerical simulations always maintain the exact cancellation of valence and ghost de- 
terminants, because rather than having the determinants of two types of quark cancel, the 
valence and ghost determinants are never introduced in the first place. Therefore, at this 
point the question arises whether we truly have a paradox. The answer is that the conflict 
between the effective and underlying theories is a real one. Given that no mixed mass terms 
ever exist in the actual mixed-action setup, we are always in the e-regime for the Goldstone 
pions of Eq. ( I3.15p . The correct prescription in this regime is to first calculate Feynman 
diagrams using chiral perturbation theory for a given orientation of the condensate, and 



8 For example, n^ s is made of a valence anti-down quark and a sea up quark, and irj v is made of a sea 
anti-up quark and a valence down quark. Not all fields are independent due to the constraint sdet(S) = 1, 
which excludes the "super-singlet" field [29| from our effective theory. 

9 The notation U(f) sca _ va i refers to a U(l) transformation with opposite phases in the sea and valence 
sectors. The diagonal U(l) sea=va i, where these phases are equal, acts trivially on meson fields. 

We checked this by explicit calculation. 
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then to integrate the result over all possible orientations. Let us return to our example of 
the two-point function (7r+(0)7r~ (x)), but now in finite volume. On the vacuum (I3.10p . its 
leading-order (LO) value will be a superposition of exponentials, as discussed above. All 
other orientations of the vacuum may be obtained by the combination of SU(2) val and U(l) va i 
rotations. These rotations leave the operators 7rf s invariant, hence their two-point function 
is unchanged when we integrate over all orientations of the mixed condensate. The finite- 
volume two-point function of pure-sea pion fields would therefore depend on m v also when 
we are in the e-regime for the Goldstone pions of Eq. (I3.15p . This prediction of the chiral 
effective theory is indeed in direct conflict with the very setup of a mixed-action numerical 
simulation. 



IV. POSSIBLE PHASES 

If the ChPT description of QCD with a mixed action is not to break down, there has 
to be some mechanism that excludes the mixed phase of the effective theory. This section 
shows that this is indeed the case: Spontaneous symmetry breaking is entirely confined to 
the sea sector. 

In Sec. IIV Al we consider the valence and ghost sectors of mixed-action QCD for arbitrary 



number of flavors. Following Vafa and Witten [24J, we employ a bound on the spectrum 
of the chiral Dirac operator to prove that none of the flavor symmetries of the valence- 
ghost sector can break spontaneously for m v ^ 0. This rules out, in particular, a mixed 
condensate. 

In Sec. IIVBI we study how this information is communicated to the effective theory. We 
begin by deriving a mass inequality in the underlying theory, which constrains the mass of 
mixed charged pions to be not smaller than the minimum of the masses of pure-sea and 
pure-valence charged pions. Specializing (for technical reasons) to the case of two flavors of 
valence quarks, we infer from the mass inequality another inequality that must be satisfied 
by the LECs of mixed-action ChPT. The LEC inequality, in turn, excludes the range of 
values that produced the paradox of the previous section. 

Finally, in Sec. IIV CI we use the results of the first two subsections to conclude that the 
only nontrivial phase structure occurs in the sea sector, where our analysis reduces to that 
of Ref. 0. 



A. Absence of flavor symmetry breaking in the valence-ghost sector 

In this subsection we prove that the full valence-ghost flavor symmetry group is not 
broken spontaneously for m v ^ 0. The analysis is carried out in the underlying theory, 
mixed-action QCD. We first consider the valence sector alone, and then extend the result 
to include the ghost sector too. 

In order to avoid cumbersome notation we consider a mixed-action theory with two chiral 
valence quarks u and d, with masses m u and m^. 11 At this point we assume that m u and 
md are both nonzero, but not necessarily equal. The valence-sector action is 

S = uDu + m u ll(l - \D)u + dDd + m d d(l - \D)d , (4.1) 



The proof generalizes trivially to any number of chiral valence quarks. 
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where D is a lattice Dirac operator satisfying the Ginsparg- Wilson relation 



30 



{ l5 ,D} = D l5 D . (4.2) 



Following Ref. 24], we consider the isosp in-breaking condensate 



{nil - \D)u - d{\ - \D)d) A (4.3) 
= ~y Tr ((1 - \D){D{\ - \m u ) + m^ 1 - [m u -> m d ]) . 

The subscript A indicates a fixed gauge-field background. Assuming that D is 75 hermitian, 
= 75D75, it follows from Eq. ( 14. 2 p that Z) is normal. Thus, .D and .D* have a simultaneous 
set of eigenfunctions with eigenvalues A and A* respectively. Once again using Eq. ( 14. 2 j) it 
follows that A + A* = A* A, so we may write 

A = 1 - e i0 . (4.4) 

Moreover, if if) is an eigenfunction with eigenvalue A, 75^ is an eigenfunction with eigenvalue 
A*; hence, the eigenvalues (|4.4j) come in pairs ±0 (except possibly at the isolated points 
= or = 7r). 

Integrating over the gauge field, 12 it is now straightforward to show that the isospin- 
breaking condensate is equal to 

- \D)u-d{\ - \D)d) (4.5) 



m 2 + 4 tan 2 |0 m 2 + 4 tan 2 |( 



where p(0) is the spectral density. 13 This clearly vanishes for m u — md — > 0, as long as 
the common value m u = irtd is nonzero. The conclusion is that there is no spontaneous 
symmetry breaking of isospin symmetry within the valence sector. 

The previous analysis easily extends to the full valence-ghost sector. Since Grassmann- 
valued condensates cannot occur (App. IA~j) . this leaves us to consider a graded-symmetry 
breaking condensate of the form 

{q v {l-\D)q v + q\l-\D)q) , (4.6) 

where q v is a valence quark and q is a ghost quark. Note the plus sign between the valence 
and ghost quark bilinears. This sign is consistent with the graded symmetries in XJ(n\n)y 
(with n the number of valence quarks), which would be broken if this condensate developed 
a nonvanishing expectation value. 14 

The condensate (14. 6p has a spectral representation similar to Eq. ( 14. 5p . with m u on the 
right-hand side replaced by the valence-quark mass m v , and rrid replaced by the ghost-quark 
mass nig, which we temporarily take to be different from m v in order to study symmetry 



12 The integration measure is nonnegative, see below. 

13 Note that p(<fi) = p{—<p). For a similar expression for the chiral condensate with overlap fermions, see 
Ref. 



31|. 



14 Recall that ghost-quark fields commute with each other. For reference, the singlet condensate that does 
not break U(n|n)y has a minus sign between the valence and ghost terms. 
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breaking. 15 Again, in the limit m g — > m v this condensate vanishes, for the same reasons as 
before, and we conclude that the full valence-ghost flavor symmetry group is not sponta- 
neously broken. 

A corollary is that no mixed condensate can ever occur. A bilinear sea-valence op- 
erator transforms in the fundamental representation of the valence-ghost symmetry group. 
Since this group does not break spontaneously, bilinear sea-valence operators cannot acquire 
nonzero expectation values. 16 

We end this subsection with a technical comment. In order to probe isospin breaking in 
the sea sector (where it can occur), we may have to turn on an (infinitesimal) difference Am 
between the masses of the up and down sea quarks. For nonzero Am, the two-flavor Wilson 
determinant is no longer positive. However, since this only happens for nonzero Am and at 
a nonzero lattice spacing, the effect is of order aAm. In this paper we work to order m ~ a 2 
only, and so we may neglect such effects. 17 



B. Absence of mixed-phase Goldstone bosons 



While in the previous subsection we showed that all is well in the underlying theory, the 
puzzle concerning the phase diagram of mixed-action ChPT remains to be resolved. We 
will begin by deriving a mass inequality relating pure-sea, pure-valence, and mixed pions in 
the underlying theory. We will then infer an inequality between the LECs of mixed-action 
ChPT, which, as promised, excludes the region of parameters that gave rise to the paradox 
of Sec. IIIIBI As in previous sections, the ChPT-level analysis is restricted to the case of 
two valence quarks. 

We begin with the following inequality in mixed-action QCD: 

tx((si d (x,y) - SUx,y)) {S sd {x,y) - S vd {x,y))) > , (4.7) 

in which Sik(x,y), i = s,v is the sea, respectively, valence quark propagator. The second 
index k = u,d denotes flavor, up or down. 18 

Although as it stands inequality (14 .7p depends only on the d propagator, the u propagator 
will be encountered shortly. In the isospin-symmetric phase, the up and down propagators 
are equal, and Eq. (" 14 . 8 [) below follows from 75-hermiticity of the Wilson operator. 

When we get back to the effective theory, we will make use of the inequality only in the 
symmetric phase. Interestingly, only little extra effort is needed to extend the inequality to 
the phase with broken isospin (the Aoki phase), so let us make this small detour. In order 
to account for the latter possibility, we add a "twisted" mass term of the form ^q s i^5T^q s 
to the (sea) Dirac operator, where q s = (u s , d s ) T . This accomplishes two things. First, it 



15 Recall that m g is necessarily positive. Assuming an even number of valence quarks we may take m v > 0, 
because, if m v < 0, we have m v — > \m v \ under a nonanomalous chiral rotation of the valence quarks. 

16 The same argument excludes a valence-ghost Grassmann condensate. However, the latter was already 
ruled out by the general proof of App. [A] 

17 For a discussion of why the argument of Ref. 24j does not apply to an isospin-breaking condensate in the 
sea sector, in which the quarks fields are of the Wilson type, see Ref. [I9I ]. 

18 While the correlation functions under study depend on two valence flavors that we have conveniently 
denoted up and down, the total number of valence quarks can be any n > 2. 
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accounts explicitly for isospin breaking. Second, it aligns isospin breaking along the third 
direction in isospin space, so that the relevant condensate, if it forms, would be proportional 
to g s 75T 3 g s . We now have that 

S ld( x i v) = ^Ssuiy, x)7 5 , (4.8) 

also when /i ^ 0. Eq. (14. 8p holds in any finite volume, and therefore also in the thermody- 
namical limit where \x is eventually turned off. 

Using that Sl d (x,y) = jsS vu (y , x)js for the chiral overlap propagator as well, inequal- 
ity ( 14. 7p can now be rewritten as 



G ss {x, y) + G m (x, y) > G sv {x, y) + G vs (x, y) , (4.9) 

where 

Gij(x,y) = (lii(x)i'Ysdj(x)dj(y)i'ysUi(y)) . (4.10) 

There are no disconnected contributions. (This is true even if we are inside the Aoki phase, 
because, with our choice of the twisted-mass term, any isospin-breaking condensate must lie 
in the r 3 direction.) 

If inequality (14. 9p holds in mixed-action QCD, it must also hold in the low-energy effec- 
tive theory, in which G ss corresponds to the sea-pion propagator, G vv to the valence-pion 
propagator, and G vs and G sv to the mixed-pion propagators. Technically, the translation 
into ChPT is done by coupling the mixed-action QCD lagrangian to pseudoscalar sources 
for the operators Ui{x)^dj{x) and dj(y)^Ui(y) [9-11, 32]. While the underlying theory is 
nonunitary, within the effective theory the decay rates of the relevant correlation functions 
are thus interpreted as pion masses. 

The LO result can be expressed in terms of the component fields of Eq. (13. 13j) as follows 

(nts( x )n7s(y)) + (^vv{ x )^vv(y)) > (^sv( x )Ks(y)) + (^s( x )^(y)) ■ ( 4 - n ) 

Here each two-point function is a tree- level (i.e., free) propagator in the effective theory, 
with mass determined by the potential, Eq. (12. 2p . For a meson of mass M, the tree- level 
propagator is 

in which r = \x — y\. Substituting this into Eq. (14. lip yields the inequality 

MssK^Mssr) + M vv K x (M vv r) > 2M sv K 1 (M sv r) , (4.13) 

where we used that M vs = M sv . If all masses are strictly positive, for large r we can use the 
asymptotic behavior of Ki, 

Kl{z) ~ {Yz e ~ Z ' N ^ °° ' (4 ' 14) 

finding 

M sv >mm{M ss ,M vv ) . (4.15) 

A nontrivial consequence of the mass inequality (I4.15P is that the LECs appearing in 
Eq. ( 12. 5p . too, are subject to an inequality. To see this, we take m s and m v large enough 
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that no spontaneous symmetry breaking takes place, and the curvatures at the origin of field 
space, which are given explicitly by Eq. (13. 8p . are all positive. We also choose m v such that 
M vv = M ss , which implies that m v = m s + 4:c' 2 . We now find that the mass inequality (I4.15P 
translates into the LECs inequality 

2 Cl -c 2 = 2c[ - c' 2 > . (4.16) 

LECs are, by definition, independent of the quark masses. While we have derived the 
inequality by considering special values of the sea and valence masses, it must therefore hold 
true for arbitrary values of m s and m v . 

We note that in our example leading up to Eq. (I3.12p . we assumed that c' 2 > 0, so that, 
with Eq. (I4.16p . c[ + 2c' 2 > c[ — 2c' 2 > 0. Therefore, c[ + 2c' 2 < never occurs, and the 
putative vacuum solution Eq. ( 13 . 1 2[) is never encountered. 



C. The phase diagram 

We now return to the ground state (Eqs. (13. ip and (I3.4p ). i.e., to the minimization of the 
potential (I2.5p . but subject to the constraints we have inferred from the underlying theory. 19 
In Sec. IIH Al we found that S g = 1. Using the result of Sec. II V A I we conclude that S OT = 1 
as well, while using the result of both Sec. IIV Al and Sec. lIVBl it follows that £ s „ = J2 VS = 0. 
In addition, we know from App. |A] that u — uJ — 0. Therefore, only S ss can take on a 
nontrivial value, with det(S<, s ) = 1, so that E ss e SU(2). Substituting the vacuum solution 



^ss 











1 











1 



S vac = 10, (4.17) 



into the potential ( 12.51) . it reduces to 

^ = -m s tr(S ss + SL)-ic 2 (tr(S ss + SL)) 2 , (4.18) 

where c 2 is defined in Eq. (13.91) . This is precisely the potential which was found in Ref. 1^ 
for QCD with two Wilson flavors, and no valence sector. Depending on the sign of c' 2 , either 
isospin and parity are broken in the sea sector (if c' 2 < and \m s \ is small enough), with 
7rf s the corresponding Goldstone bosons, or there is a first-order phase transition (if c' 2 > 0), 
with a minimum nonvanishing pion mass at nonzero lattice spacing. 

In summary, adding a chiral-fermion valence sector to the dynamical Wilson-fermion 
theory enlarges the flavor symmetry group, and, in principle, allows for many new symmetry- 
breaking condensates. Nevertheless, none of these condensates actually develops, and we 
have recovered the usual two-flavor phase diagram for the sea sector. 



19 In this paper we analyze the potential for the case of two valence quarks, but since the proof of Sec. IIV Al 
is valid for arbitrary number of valence quarks, we expect that the analysis of the potential can be 
generalized accordingly. 
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V. CONCLUDING REMARKS 



The chiral lagrangian for lattice QCD with a mixed action, with Wilson sea quarks and 
chiral valence quarks, is rather involved. We have studied this chiral lagrangian in the case 
of two sea and two valence quarks in the LCE regime, m ~ a 2 . With no restrictions on the 
values of the order-a 2 LECs, the phase structure is rather intricate. 

In particular, the effective potential appears to allow for a mixed phase with a condensate 
pairing valence with sea quarks. Such a mixed phase would contradict the very setup 
of mixed-action simulations, because pure-sea correlation functions would depend on the 
parameters of the valence sector. 

As we have shown in this paper, the underlying theory excludes such a phase. By ex- 
tending the well-known argument of Ref . [24] , which relies on the existence of a gap in the 
spectrum of the chiral Dirac operator for nonzero mass, we have shown that none of the 
flavor symmetries of the valence and ghost sectors can be broken spontaneously, for any 
number of valence quarks. This forbids, in particular, a mixed condensate. 

In addition, we have shown that the mixed-pion mass cannot be smaller than the min- 
imum of the pure-sea and pure-valence pion masses. This mass inequality translates into 
an inequality between LECs appearing in the chiral lagrangian. The inequality, Eq. (I4.16p . 
must hold independent of the choice of action of the underlying lattice theory. In terms of 
the original LECs in Eq. ( I2.2p . the bound reads 

2W M -Ws>0. (5.1) 

Without this restriction on the LECs, the effective potential by itself would allow for a 
mixed phase. The inequality we found follows from expanding the effective potential around 
the trivial vacuum, i.e., Eq. (I3.8p . and imposing the mass inequality (I4.15p . It is possible 
that more constraints on these LECs exist, that would follow from a complete study of the 
effective potential for arbitrary S vac , by imposing the constraints we found in the underlying 
theory. 



Not surprisingly, we recover the well-known conclusion of Ref. [19[ about the possible 
phase structure in the sea sector. If c' 2 > 0, there is a first order transition when m s changes 
sign, and the pion mass is always larger than zero; if c' 2 < 0, a second order transition occurs 



for small enough m s , in which isospin and parity are spontaneously broken 23]. This raises 
the question whether an inequality might also be derived for c' 2 by considering the charged 
and neutral pion masses in the sea sector. The reason that this is not possible, however, 
is that the neutral pion propagator in QCD with broken isospin contains "disconnected" 
diagrams, so that the arguments of Sec. IIVBI do not apply. Note that in Sec. HVBI we only 
considered propagators for charged pions, making sure that no disconnected contributions 

90 

appear. 

Finally, we expect that our conclusions generalize to other mixed actions with a chiral 
valence sector. For instance, if the sea quarks are staggered, they might exhibit a nontrivial 



phase structure [33(. In a mixed-action theory with a staggered sea sector, this nontrivial 
phase structure would remain confined to the sea sector, just as in the case we considered 
in this paper. 



20 Here "disconnected" diagrams are diagrams with disconnected quark loops; they are still connected if also 
gluon lines are taken into consideration. 
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Appendix A: Absence of Grassmann-valued condensates 

In this appendix we prove that Grassmann-valued condensates cannot arise. To be con- 
crete, we consider a possible (q s q v ), where q s is a sea quark and q v a ghost quark. The 
relevant part of the path integral is of the form 21 

N 

z ig) = Yl d( fi d( ii < ^i d( ii ex P [-qXq - - ^ B( i - q ] Yq] , (Ai) 
J i=i 

where X corresponds to the sea-quark Dirac operator and Y to the ghost-quark Dirac oper- 
ator. Here the variables q^ q i7 i = 1, . . . , N are independent Grassmann-valued (fermionic) 
variables, and the variables q^ are bosonic, c- number variables. The matrices X and Y 
contain c-number entries, while the matrices A and B are Grassmann-valued. Of course, 
the matrices A and B vanish in mixed-action QCD. But, in order to study the possible 
occurrence of a mixed sea-ghost condensate, one chooses (appropriate entries of) A and B 
nonzero, and then one takes the limit of A and B to zero after the volume has been taken 
to infinity. 

Let us first recall the standard case for which also the qi are fermionic. 22 An example 
is the formation of an Aoki condensate for two flavors of Wilson fermions. In this case 
we would take X = Y equal to the hermitian Wilson-Dirac operator, and we may take 
Aij = rriAdij, B^ = —m B 6ij. 23 Without loss of generality, we may assume that X has been 
diagonalized, 

We then find that 



Xij — Xi5ij . (A2) 
N a N 

£te> (9) - ~^ {q) = ~Z [q) E ~TT~ — 2 . (A3) 
^— ' am a xf + m z 

i=i i=i 1 

where in the last step we took = m# = m, and where (. ..)( ? ) is the unnormalized 
expectation value with partition function Z^ q \ After averaging over the gauge field and 
taking the infinite-volume limit, followed by the limit m — > 0, one finds a nonvanishing 



condensate if and only if there is a nonzero density of near- zero modes [34|, because 



in 

lim — = tt5(x) . (A4) 



21 The argument for the vanishing of (q v q v ), with q v a valence quark, is similar, except that in this case 
X = Y. 

22 In that case, the q\, like the q i , are independent fermionic variables as well. 

23 For vtla = ms this corresponds to an Aoki condensate pointing in the <J2 direction in isospin space. 
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Now let us return to the case in which the q~i are bosonic. The entries of A and B are 
now fermionic, and we can, in fact, use this to work out what happens in any basis. Taking 
Aij = a5a 5jj for some fixed values of %q and jo, the integral (lAlj) now evaluates to 

Z {q) = sdet ( * y J = det (X - AY^B) /det(F) (A5) 

= det(X)exp(-ti(X- 1 AY- 1 B)) /det(F) 
= det(X) (1 - triX- 1 AY^B)) /det(F) , 
and we find (taking Grassmann derivatives to be left derivatives) 

— = - {Y-'BX^). oio det(X)/det(y) . (A6) 

It is clear that in this case, in contrast to the standard case reviewed above, no nonvanishing 
value can occur for a Grassmann-valued condensate in the limit that B — > 0. The basic 
reason is that we can always expand the condensate in terms of the components of A and B, 
resulting in a finite polynomial in those components. When we take A and B to zero, the 
corresponding Grassmann-valued condensates will thus always vanish. We conclude that 
when we consider the vacuum expectation value for the field E in Eq. (12. ip . we can set 
lo = uJ = 0. 

Two more comments are appropriate. First, one might wonder what would happen if one 
takes the matrices A and B to be bosonic. In this case, Q = qAq + q'Bq is Grassmann- 
valued, and it is easy to see, by expanding the exponent in Eq. flAip in terms of Q, that in 
this case the partition function does not depend on A and B at all. Condensates would thus 
trivially vanish. 

Our second comment is that the same question can also be studied in ChPT. In other 
words, one can assume that a priori u and uJ in Eq. (12. ip do not vanish. One then finds 
that the equations of motion (for constant fields) in the effective theory dictate that oj and 
W vanish. This is, of course, consistent with the QCD-based argument given above. 



Appendix B: Proof of Eq. Q3T5j) 

Consider the matrix (c/. Eq. OIL 



£ q = ( .5," ^ sv ), (Bl) 

V ll/ - I vs ^vv J 

where all the entries on the right-hand side are two-by-two matrices with complex entries. 
As we have seen, E q is an element of SU(4). The unitarity constraint E q £t = 1 provides 
constraints on the submatrices: 

£ ss £t s + x gv xl = i ; (B 2a) 

E vs El + S w EL = 1 , (B2b) 
S SU E OT = S SS S^ S . (B2c) 

The effective theory is invariant under independent flavor rotations in the sea and the valence 
sector. With V s € U(2) sea and V v G U(2) va i, S q transforms into 

( V S -L 8S V} iV^ s X \ m ^ 

^ \iv v z vs vj KS w K f J ' 1 } 
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We now want to prove that this matrix can be brought into the form (I3.5P if E OT = p v l, 
with p v an arbitrary complex number. 
First, from Eq. (1B2bj) we conclude that 

Ks = V^-\Pv\ 2 C, (B4) 

where C G U(2). We also find the bound \p v \ < 1. Equation (IB2cl) implies that 

SsuS OT E OT E sll = Yj ss Tj vs Tj vs Yj ss , hence 

\p v \ 2 E av El = {l-\p v \ 2 )E sa El • (B5) 
Substituting this into Eq. (lB2al) . we obtain S ss Ej s = E W E^ = \p v \ 2 , and we may thus write 

CIS 

Zss = P* v A, (B6) 



where A G U(2). From Eq. (|B5]) we obtain E sv = a/1 - \p v \ 2 B , with, from Eq. (lB2c|) . 
£> = AC^. Performing a flavor transformation with V v = C\ we can thus write E q as 

^ = ( f:t %p :t ) > (B7) 



.W 1 P^ 1 

where we introduced = a/1 — |p„| 2 . 

The matrix A can be written as exp(ia At r M ), where r M = (1,t). Hence we can define 
U = exp(—ia fl T^/2), with AU 2 = 1. Performing a flavor rotation with V s = U and writing 
C/t = e **y W ith V G SU(2) simplifies E q to 

/- p ; e ^ 2 ip^V 



E q = e * ^ ' ^« ' . (B8) 

q y zp^ p„e 4<?, y v 7 

Finally, the SU(2) matrix can be diagonalized by a flavor transformation V s = V v such that 

V -»■ D = exp(z#r 3 /2) . (B9) 

Absorbing the phase e~ % ^ into p v , we arrive at 

S q = e» ( ^ if * D ) . (BIO) 

This is precisely Eq. f 13 . 5 j) . Because E q G SU(4), it follows that det(E q ) = e 4 ^ = 1, and thus 
that = mod 7r/2. 
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